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I briefly introduce the methods by which lattice QCD predictions for 
RHIC are obtained. Next I deal with lattice determinations of strangeness 
. production and event-to-event fluctuations of conserved quantities. I also 

^SJ ' present a new diagrammatic method for computing derivatives with re- 

spect to chemical potentials, and conclude with discussions of some tests 
' of thermal perturbation theory which follow. 
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p. 



Heavy-ion collisions at the RHIC have already given evidence for dense 
and hot matter |||], and may lead to a discovery of the predicted plasma 
^ ' phase of QCD if such predictions are made precise enough. In recent years 

computations in lattice field theories have become precise enough to con- 
front phenomenological analyses of experimental results. Among the many 
interesting results from RHIC Q I single out three for comment — indica- 
tions of early thermalisation leading possibly to hydrodynamic flow, rapid 
chemical saturation of strangeness, and fluctuations from one event to an- 
other. Each of these observations may be related to quantities which are 
easily computed in finite temperature lattice QCD. Flow is strongly con- 
nected to the equation of state, strangeness to the Wroblewski parameter, 
and fluctuations to various susceptibilities. The equation of state has been 
adequately dealt with elsewhere |p, and I will restrict myself to the rest. 

Interesting thermodynamical quantities can be constructed by taking the 
derivative of the free energy with respect to intensive quantities. Consider 
QCD with the intensive variables temperature, T, and the quark chemical 
potentials, fiu, jJ-d-, IJ-s (for flavours u, d and s of quarks). The first derivative 
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of the free energy, F, with respect to one of the chemical potentials is the 
quark number — 

d d 

{nf) = ^ — F{T,^iu,fi'd,fi's) = ^ — log Z{T, flu, ^J'd:^J's), (1) 
dfj,f dfif 

where Z is the partition function. The second derivatives are called quark 
number susceptibilities Q 

In general such second derivatives measure microscopic fluctuations in equi- 
librium. It has recently been discovered that these fluctuations, Xfg-, may 
be directly accessible in heavy-ion collisions j^. Higher derivatives, which 
we deal with later, may be called non-linear quark number susceptibilities, 
in analogy with similar quantities in condensed matter physics. 
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Fig. 1. Finite temperature Euclidean field theory in the continuum and on a lattice. 

I will describe information on these derivatives which we have obtained 
from lattice simulations of QCD. These are numerical estimates of the QCD 
partition function, Z, in the Euclidean thermal field theory context — 

Z(T,Hu,tJ-d,tJ's) = J dU Yl detM{mf,Hf)exp[-S{U)]. (3) 

f=u,d,s 

In this formula S{U) is the gauge part of the action, M is the Dirac operator 
for quarks of mass nif with chemical potential fif, and the integration is 
performed over all configurations of gauge fields. The chemical potentials 
enter the Dirac operator as if they were constant U{1) gauge fields (see 
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Figure |^). In the Euclidean formulation of finite temperature field theory 
the temperature enters indirectly through the fact that the Euclidean 'time' 
direction is of extent 1 /T Q . 

Since such integrals have the usual ultraviolet divergences of field theory, 
they can be defined with a space-time lattice as a regulator. The spacing 
between lattice sites in all directions is a. The gauge field associated 
with infinitesimal changes in position, i.e., d^, is replaced by the finite 
transporter exp(m^^). The number of lattice sites in the time direction is 
Ni. At fixed temperature this gives the relation 

aNt = 1/T, (4) 

This is used to eliminate the ultraviolet cutoff, A = l/o from all computa- 
tions in favour of the physical scale T. After this is done, the regulator must 
be removed by taking A ^ oo while holding fixed all physical quantities. 
This process is called "taking the continuum limit", since it means that 
a —5- at fixed temperature by taking Nt oo while holding T constant 
(see Figure ||) . 

Taking the continuum limit is exactly the same as normalizing the field 
theory. A measurement of the renormalized strong coupling at the scale of 
1/a flows according to the two-loop /9-function of QCD. As a result, good 
control over the continuum extrapolation comes from perturbation theory, 
yielding [§ 

^^ri. 15 ±0.05 {Nf = 0), 

A^ ~ 1 0.49 ± 0.05 {Nf = 2). 

This gives us good precision in pinning down the running coupling at any 
given temperature since log{T/K-^) = \og{T/Tc) + log(rc/A^). Tc/Aj^ 
is a reasonably easy quantity to measure because there are nice definitions 
of the renormalized QCD coupling on the lattice which show the usual log- 
arithmic scaling without any power corrections in a. Other quantities may 
have power corrections which need to be subtracted before the logarithmic 
scaling can be seen. While this is tedious, the great advantage of the lattice 
is that it allows full control over infrared divergences which plague finite 
temperature field theory. 

Another fact is crucial. In a lattice computation we do not determine 
the integral in eq. (^) before taking its derivatives. Instead, we take the 
derivatives before doing the integral numerically. For example, we notice 
that for any matrix M{x) where each matrix element may depend on some 
variable x 

^^^^ = detM(x)TVM'M-\ (6) 
ox 
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where M' denotes the matrix each term of which is the derivative of the 
corresponding term of M. As a result, 



(n. 



) = 1/ dUTTMiM-^\{MfeM-S{U)] = (TVM/M-i). (7) 



The expectation value on the right is computed with a Monte Carlo proce- 
dure which simulates the integrand of eq. (^) [^. Such a procedure works 
only when the integrand is non-negative. Since the Euclidean Dirac opera- 
tor with chemical potential has complex eigenvalues, the determinant is not 
positive definite, and lattice Monte Carlo simulations of QCD at finite chem- 
ical potential become tremendously hard to do. In all the work reported 
here we deal with the susceptibilities evaluated at zero chemical potential. 
Interestingly, they can be (and have been) used to continue lattice QCD 
information to non-zero chemical potential Q. 

Finally a word about flavour symmetry breaking. If the u and d quark 
masses in nature were equal then flavour symmetry would be broken only 
in electro-weak interactions. Lattice QCD computations are done in this 
limit. In reality, however, u and d quark masses differ. It turns out that 



this breaking is almost irrelevant to thermodynamics |1C]. The strange 
quark is much heavier, with a mass not much different from Tc- Hence it is 
almost quenched close to Tc but should be treated as unquenched far above 



Tc. 



2. Fluctuations 

2.1. Lattice Measurements 

I will introduce some notation. The usual baryon chemical potential is 
/^o = (/"m + /"d + /"s)/3. The chemical potential associated with the isospin 
quantum number is = (/i^ — fj-d)/"^- The corresponding number densities 
are (no) = {uu + rid + '^s)/3 and (723) = (n„ — n^)/2. These are zero 
whenever the chemical potentials vanish. The susceptibilities obtained by 
taking double derivatives with respect to fiQ and /U3 are written xo and xs- 
These can be non-zero even when the chemical potentials vanish. 

We recently improved upon previous measurements |ll| of these quan- 
tities in several ways. First, by changing the size of the spatial box within 
which the lattice computation is done, we have found a range of sizes such 
that the box has no effect on the physical measurement. All our subsequent 
measurements are in this range of box sizes. Secondly, unlike previous com- 
putations, we have held the quark mass, m, fixed in terms of physical mass 
scales such as GeV or Tc as we change the temperature. Previous studies 
had, for convenience, fixed ma = m/TNt, which meant that their quark 
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Fig. 2. The susceptibility X'i s-s a function of T for several different valence quark 
masses, m.y and lattice spacing a ~ 1/4T. The lines summarize data for quenched 
QCD while the symbols are for data in Nf = 2 QCD for the m^/Tc values shown. 
Error bars are mostly smaller than the symbols. The shaded region covers a range 
of quark masses appropriate to strange quarks. 



mass changed as they scanned across temperature. Finally, by improving 
the estimators of the traces, it turns out that modern computers can allovi^ 
us to reduce the error bars in some of the susceptibilities by over 3 orders of 
magnitude — an effective gain of a factor of a over milhon in the statistics 
available ten years ago when the last computations were performed. 

For a = l/4r, i.e., Nt = 4, our simulations with the quenched theory 
[Nf = 0) ||l^ and with two light flavours of sea quarks {Nf = 2) |13| showed 
that the susceptibility is quite different from that for the ideal gas, Xfft, 
on the same lattice. There is a small effect from unquenching the light sea 
quarks — about 5%. Since this effect is so small, it seems that the effect of 
unquenching the strange quark should be smaller than the statistical errors 
in our measurement. By making the valence quark heavier we can therefore 
investigate the dynamics of strange quarks. We will return to this important 
point later. 

The quark mass dependence of xs is quite nontrivial and is interesting 
in itself. It has been known for a long time that X3 can be written as 
the zero-momentum limit of a certain component of a vector correlation 
function at finite temperature @ . Now, the breaking of Lorentz symmetry 
at finite temperature due to the selection of a preferred frame (that of the 
heat bath) means that angular momenta do not necessarily label the states 
of the system [^]. One of the most well-known consequences of this is 
the difference between electric and magnetic polarizations of gauge bosons 
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Fig. 3. The dependence of X3 s^nd the pseudo-scalar susceptibihty, x^r, on the 
common screening mass of the scalar and pseudo-scalar at T = l.STc and 2Tc. 
Error bars are smaller than the symbols. Xrr is the zero- momentum pseudo-scalar 
correlator lIT3. 



at finite temperature [jlGf. In any case, this broken symmetry causes xs to 
mix with scalar/pseudo-scalar representations [0, |l^. Previous work on the 
lattice has seen clearly that correlations in this channel at high temperatures 
deviate strongly from that of an ideal quark gas ||T^. As the quark mass is 
changed, this correlation function, the corresponding screening mass, and 
X3 all change in response (see Figure |3|) . 

The flavour off-diagonal susceptibility Xud turns out to be surprisingly 
small. Our measurements reveal that above Tc the dimensionless number 
Xud/T'^ is zero to within a few parts in 10^. This is a major surprise, because 
a recent computation in resummed finite temperature QCD shows that this 
quantity should be of the order of a|(27rT) log as(27rT) and predicts that 
it should be of order 10~^ [19|. A non-log contribution of order a^{27rT) 
remains to be computed, but even if this cancels the computed term at some 
r, the range of temperatures over which results are available is large enough 
that a substantial non-zero value would still be seen. This disagreement 
between the lattice and perturbative computations stand as a puzzle. 

Below Tc this off-diagonal susceptibility has only been investigated in 
the quenched theory. It is small but clearly non-zero (see Figure ^). With 
changing quark mass it is seen to vary roughly as where m^r is the 

pion screening mass, showing that such fluctuations are essentially carried 
by pions. The connection between the vector- vector correlator Xud and the 
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Fig. 4. The dependence of Xud on TO,r at T = O.TSTc in quenched QCD. 



pion below Tc comes from the fact that the contribution 



(8) 



to Xud is allowed at finite temperature, and therefore, dominates p exchange 
contributions purely kinematically. This is another manifestation of the 
same physics that led to the correlation shown in Figure |^. While Xud is 
non vanishing below Tc, X3 is consistent with zero. 



Continuum extrapolation of these results have been attempted |20, 21 



22 1 . It turns out that X3 measured with staggered Fermions have large 
power corrections in a. As a result, the ratio Xs/Xfft at T = 2Tc decreases 
from a = 1/4T to 1/6T but then turns over and approaches the limit from 
below. The same limit is obtained by extrapolating Xs/^^ using the usual 
staggered Fermions or an improved version, although both these extrapola- 
tions approach the limit more smoothly. The results at 2Tc are in marginal 
disagreement with the resummed perturbative computations of [|19[ (i.e., 
disagrees at the 1-cr level but agrees at 3-0"), but become compatible with 
it at r = 3Tc. Xud/T'^ remains compatible with zero at the level of a few 
parts in 10~^ in the continuum limit. 



2.2. Applications to phenomenology 

Two important qualitative observations emerge from the lattice com- 
putation. First, that above Tc one has non-vanishing X3 but Xud is zero. 
Second, that below Tc xs vanishes and Xud is non-zero. Since there are only 
two independent types of susceptibilities (as I show later), all fluctuations 
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of interest are governed by these two and their changes with valence quark 
mass. 

Fluctuations of electric charge, for example, are controlled by the sus- 
ceptibility 

X<7 = g (10X3 +Xs + Xud - 2Xns) , (9) 
whereas baryon number fluctuations are related to 

Xo = ^ (4X3 + Xs + ^Xud + 4x«s) . (10) 

The following quantitative conclusions can be obtained — 

1. For T » Tc, since X3 ~ Xs » Xud ~ Xus, we have Xq ~ (11/9)X3 and 
Xo ~ (5/9)x3 so that the ratio Xq/XQ ~ 2. 

2. When T > but very close to Tc, since X3 > Xs » Xud ~ Xus, we 
find Xq ~ (10/9)x3 and xo ~ (4/9)x3 so that the ratio Xq/Xo ~ 2.5. 

3. For T < Tc since Xud l/rn^ and assuming that Xus oc l/rnj^, since 
X3 ~ Xs ~ 0, we expect that Xq/Xo ~ 0.25 + 0{ml/mj^). 

Under the assumptions given above, there are diametrically opposite 
predictions above and below Tc — 

XO < Xg < Xs (T > Tc) 

XO>Xq> Xs {T < Tc) (11) 

The ordering of fluctuations in baryon number |^^, charge and total 
strangeness are therefore radically different above and below the phase 
transition. 



3. Strangeness production 

Strangeness abundances in heavy-ion collisions at the CERN SPS col- 
lider and the RHIC have been analyzed extensively. There is some consensus 
that the observed chemical composition is that in equilibrium close to Tc, 



and that it cannot arise due to hadronic rescatterings |25]. One of the central 



quantities that has been extracted from data is the Wroblewski parameter 

(nn) + {da) 



The averages on the right are defined to be the number of primary created 
quark pairs of each flavour. It turns out that lattice determinations of static 
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equilibrium quantities can be used to predict this dynamical quantity under 
some well-defined, and testable, assumptions. 

As a preparatory example, consider the electrons in a metal interacting 
with external fields. In a static magnetic field, if, at a fixed temperature, 
the response is an induced magnetization whose rate of change with the 
field strength is the magnetic susceptibility, x(0)- This is a measure of the 
fluctuations of spins in thermal equilibrium. On the other hand, when an 
electromagnetic wave of frequency a; propagates through the medium, it is 
attenuated due to dissipative phenomena which generate many excitations 
in the medium. One can describe the dissipation through a complex sus- 
ceptibility x('-^)) describing the response of the material to a magnetic field, 
H{uj) of frequency u> ||2^. Causality relates the real and imaginary parts of 
xi^) through a Kramers-Kronig dispersion relation. From the fluctuation- 
dissipation theorem it is possible to deduce that the complex susceptibilities 
are proportional to the static susceptibility if the characteristic time scales 
of the system are very different from the energy scales dominating the pro- 



duction process |28|| . 



Quenched QCD (T^,) 
RHIC Au-Au 
SpS S-S 
SpS S-Ag 
SpS Pb-Pb 
AGS Au-Au 



AGS Si-Au . 9 

0.^2 0^3 ijtB 0^ 



Fig. 5. The Wroblewski parameter on the lattice compared with extraction from 
data The error bars are statistical errors. For the lattice extraction the 

bracketed interval an estimate of possible errors due to extrapolation to Tc. 

This carries over to strangeness production. The rate of production of 
quark pairs is proportional to a complex susceptibility, and hence to the 
static susceptibility that we measure. This gives 

A. = = ^, (13) 

Xu + Xd Xu 



where the susceptibilities are evaluated at the temperature and chemical 
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potential characteristic of the cohision. In Figure ^ we display our prediction 
of Xs, from the lattice computations already outlined, and a comparison with 
values extracted from experiments. Further details, including a complete list 
of all assumptions and ways to test them, can be found in [pO|. 



4. Testing perturbation theory 

We have already seen evidence of non-perturbative effects in the suscep- 
tibilities. We extend this observation by investigating non- linear suscepti- 
bilities. We define these as higher derivatives of the free energy — 

Q"log^ .... 

where the flavour indices /, g, h, etc, need not be distinct. There is a nice 
and systematic way of evaluating the derivatives of Z. It begins by noting 
that the chemical potentials Hf appear in the partition function of eq. (|3|) 
only through the quark determinant. Then, we can evaluate the derivatives 
at by a chain rule starting with 



ddetM , , _ (1) a^detM 
/ 



„ det MOy>, -——— = detM 



0f^0^^^+6f,0f], (15) 



etc. The operators Oj are defined recursively through the relations 

fl/-o('^-i) 

of ^ 

and the concrete computational rules are — 



(U , 1 OMf , , 

Oy =TiM'jMj^, ^ = -Mj^M'fMj^. (17) 



This is the complete set of rules for writing down the operator expressions 
for the non- linear susceptibilities [ |29| . 

In the continuum theory, since the Dirac operator contains /x linearly, 
second and higher derivatives, M'j etc., vanish. Every Afj- corresponds to 
an insertion of 70 (where Aj is a flavour generator) |Q, and each M^^ 
is a quark propagator. Thus the chain rule (eqs. [T5|-p^) can be written 
diagrammatically. The rules for a susceptibility of order n are — 

1. Put down n blobs (each corresponding to an Mj, i.e., a derivative 
with respect to ^u/) and label each with its flavour index. 
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Fig. 6. The operator topologies which contribute to the third and fourth order 
susceptibihties. 



2. Join the blobs by lines (each representing an MJ ) into sets of closed 
loops such that each loop contains only blobs of a single flavour. Count 
the number of ways in which each topology arises and sum them all 
up. 

3. For degenerate flavours ioi fif = 0, the operator depends only on the 

(n) 

topology and the flavour label on Ojr is irrelevant. So delete all the 
flavour indices after the counting is done. 

4. The operators can then be labeled only by the topology, which is 
specified completely by the number of blobs per loop and the number 
of such loops. Thus, each distinct topology is a partition of n. 

In Figure ^ are shown the topologies that contribute to the 3rd and 
4th order susceptibilities. There are clearly two types of operators — one 
quark line connected operator for each n, and the remaining quark line 
disconnected. Flavour off-diagonal operators are necessarily quark-line dis- 
connected. Since the free energy is even in each Hf, the odd-order suscep- 
tibilities vanish for /i = 0, just as do the number densities. 

The rules show that Vxuu/T = {O2 + di) and Vxud/T = (On). Due 
to flavour symmetry, Xuu = Xdd- The number of independent physical 
quantities, i.e., susceptibilities, is equal to the number of operators. Hence, 
the operator expectation values are themselves physical. At third order 
Vxuud/T = {Oui + C12) and Vxuuu/T = (Om + 3012 + O3). Flavour 
symmetry gives two different physical quantities — Xuuu = Xddd and Xuud = 
Xudd whereas there are three different operators. At fourth order, there 
are three different physical quantities, which are Xuuuu = Xdddd, Xuuud = 
Xuddd and Xuuddi but flve different operators. Due to flavour symmetry, the 
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number of different n-th order susceptibilities is equal to 1 + n/2 for even 
n and (1 + n)/2 for odd n, whereas the number of distinct matrix elements 
is the number of partitions of n. For n > 2 there are more matrix elements 
than susceptibilities, and the former cannot all be physical. The particular 
fourth order susceptibility — 

Xuudd = ^(Ciiii + 20211 + O22) - xlu - '^Xld (18) 

is a cumulant related to the Binder variable and hence interesting to 
study. 

For each n only one of the susceptibilities, that with only a single flavour, 
contains a quark-line connected diagram. All other susceptibilities are nec- 
essarily quark-line disconnected. We have investigated some of these quark- 
line disconnected quantities numerically. In dynamical QCD with Nf = 2 
at temperatures T > 1.5Tc, it turns out that Xud/T"^ vanishes to one part 
in 10^, and both Xuud/T'^ and Xuudd/T^ vanish to better than one part in 

While the quark-line disconnected diagrams are expected to vanish in 
an ideal gas, in QCD they may be connected by gluon lines, and dressed 
by all possible gluon and quark loops. In certain power counting rules 
were developed which may be applied to operators such as these: the main 
ingredient being that every loop with n blobs connects to Ug electric gluon 
lines, where Ug > 1 and -|- n is even. As a result, (Cii) oc (actually 
In g as shown in |jl^ after a detailed computation) . All contributions to 
the third order susceptibility vanish. Of the diagrams contributing to Xuudd, 
{O22) oc (7^ and gives the leading perturbative contribution. At tempera- 
tures of 2Tc, for Nf = 2, we get {Ou)/T'^ ^ 0.1, and Xuudd/T^ ~ 0.5. These 
rough perturbative estimates can easily be modified by an order of magni- 
tude due to subleading logarithms and numerical coefficients. Nevertheless, 
the lattice results are significantly below the perturbative estimates, and 
temperature independent over a range of temperatures where the perturba- 
tive estimates vary by a factor of 5. 

This finite temperature analogue of Zweig's rule holds in a region of 
temperatures away from Tc- Closer to there there is some evidence for 
non-zero values of Xud H, 11, ^] as well as Xuudd- Since these quantities 



measure departures from ideal gas behaviour, they would be very interesting 
quantities to study in the vicinity of the critical point of QCD. 

It is a pleasure to thank the organizers for a wonderful school. I would 
also like to thank Jean-Paul Blaizot and my collaborators, Rajiv Gavai and 
Pushan Majumdar, for discussions. 
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